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An approximation method to study the properties of a small black hole located on the TeV brane 
in the Randall-Sundrum type I scenario is presented. The method enables us to find the form of the 
metric close to the matter distribution when its asymptotic form is given. The short range solution 
is found as an expansion in the ratio between the Schwarzschild radius of the black hole and the 
curvature length of the bulk. Long range properties are introduced using the linearized gravity 
solution as an asymptotic boundary condition. The solution is found up to first order. It is valid in 
the region close to the horizon but is not valid on the horizon. The regularity of the horizon is still 
under study. 



o 

O ' I. INTRODUCTION 
(N 

, Within the Randall-Sundrum brane world scenarios P, |2| the four dimensional universe is viewed as a three-brane 

■ with brane tension cr, embedded in a AdS five-dimensional bulk, with a cosmological constant A. A flat brane is 
[ achieved by fine-tuning the brane tension and the cosmological constant 
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SttGsA^-^ ; STrGsa-i- , (1) 

m ' 

^ . where £ is the curvature length of the AdS bulk. 
\^ ' The Randall-Sundrum I (RSI) scenario was postulated Q to solve the hierarchy problem by making gravity strong 
at the weak scale (of a TeV). In this scenario the world consists of two opposite tension branes embedded in an 
ADS bulk with an orbifold structure. The negative tension brane is called the TeV brane, while the positive tension 
\ brane is called the Planck brane. Z2 orbifold symmetry is assumed about each brane. The ratio between the TeV 
' electro-weak scale and the Planck scale is taken to be A '--^ 10~^^. In this scenario the standard model particles are 
I assumed to propagate only on the TeV brane. Gravity, on the other hand, is five dimensional and propagates in 
— «. ' the bulk. The five dimensional gravitational coupling is of order ITeV, G5 = {M^)~^, M5 ^ ITeV. The hierarchy 
appears on the TeV brane at distances much larger than the curvature length £. The gravitational attraction becomes 
I ordinary four-dimensional attraction with Newton's constant G4 ~ (Afpj)^^, Afpj = lO^^TeV. 

RSII consists of a single positive tension brane with Z2 symmetry about the brane. This model does not yield a 
. . , low scale gravitational theory. From now on we restrict our attention to the RSI scenario. 

^ ' Black holes in theories with extra dimensions have been studied widely. Myers and Perry Q found Schwarzschild 
k>( ] type solutions (MPS) in D-dimensional asymptotically flat space. Black hole solutions were also found in asymptoti- 
''^ ' cally AdS space 0, No non-trivial closed form black hole solutions, other than the black string solution Q, which 

■ extends in a uniform manner from the brane into the extra dimension, have been found in three-brane theories of 
the Randall-Sundrum type. Given that there is considerable interest surrounding the production of black holes at 
accelerators , and in collisions of cosmic rays Q , it is important to develop approximate methods to find black hole 
solutions in Randall-Sundrum brane world theories. 

Some initial attempts at finding black hole solutions centered on deriving the induced metric on the brane by 
solving the Hamiltonian constraint conditions |^ . Some of the induced solutions do not arise from matter distributions 
confined to the brane ^3 • Linearized solutions about RS backgrounds ^Tl\ as well as numerical solutions have also 
been derived. In a recent paper, Casadio and Mazzacurati .13] investigated how solutions of the Einstein equation 
in RSII models propagate into the bulk by using a asymptotic expansion of metric coefficients in r~^, where r is the 
radial distance on the brane. They reduced the Einstein equations to a set of differential equations involving the 
expansion coefficients which are functions of the bulk coordinate w. Solving these equations they found a family of 
solutions parameterized by two physical parameters, the mass M and the post-Newtonian parameter rj. 

It would be interesting and important to find black hole solutions in RSI models with TeV scale gravity. We are 
especially interested in how such solutions evolve when the black hole mass increases starting from around a TeV 
(where the flat space MPS limit is valid) to macroscopic scales. From the point of view of cosmology it is useful to 
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understand how TeV scale primordial black holes accrete or radiate matter when immersed in a hot plasma. For that 
we need to work out their thermodynamic properties. That requires a detailed determination of the area of the event 
horizon as a function of mass. As the size of the black hole approaches the AdS radius, which is of the same order 
as the length of separation of the two branes , the character of the solution should change and Gregory-Laflamme 
type instabilities [lj| could materialize. One has to find bulk profiles of TeV mass black holes to address these issues. 

A TeV black hole solution in RSI scenario is characterized by two scales: The curvature length of the bulk 
generated by the bulk cosmological constant and the five dimensional Schwarzschild radius, /i, related to the mass of 
the black hole as /i = ^/SG^M/JStt) . For a reasonable setup 15] the curvature length in RSI is of order £ ~ 10(TeV)~^, 
while the "mass" can range from a TeV black hole /i ~ l(TeV)~^, to ordinary black holes ^ ^ 10^^(TeV)~^. Clearly, 
one should use an appropriate approximation method for various values of the parameters £, fi. 

Linearized gravity is valid when the scales under study are much larger than the Schwarzschild radius r ^ fi. 
Einstein equations are expanded to first order in /i^ ~ G^M , and the solution appears to be a superposition of various 
modes; the zero mode, the massive modes, and the radion. The radion appears to be a four dimensional massless 
mode which describes the relative motion of the branes. Being massless, the radion mediates a long range strong 
gravitational force and therefore ruins the hierarchy. Provided a stabilization mechanism which in general makes 
the radion massive, the linearized solution demonstrates the transition between five dimensional TeV scaled gravity, 
at distances smaller than the curvature length, and four dimensional Planck scale gravity at distances larger than 
the curvature length along the TeV brane. The advantage of the asymptotically conformally flat linearized solution 
is to enable the definition of conserved momenta and identification of mass (in analog to the ADM mass). On the 
other hand, one is unable to study the horizon of the black hole using linearized gravity, since the horizon is located 
approximately dX r '-^ ^. Linearized gravity in RSI brane scenario was widely studied [l|, |3, ^M, ILTI lla | and will be 
presented in appendix 1X1 Asymptotic expansion in |l3j | is valid only when r <C /i. In cases where fj, < i the 
expansion can not be used for r < i, and will not enable us to study the horizon. 

In this paper we use the following method: we expand the solution in the dimensionless parameter e = fj,/£ <^ 1. 
The coordinates are rescaled by /z, therefore, an expansion in e actually means an expansion in inverse powers of £. 
Zcroth order in e should be nothing but the five dimensional flat background MPS black hole First order in e 
"turns on" the brane tension, while second order in e "turns on" the bulk cosmological constant. 

Up to first order in e (e solution) the metric is found in section ^ The e solution is required to be regular 
everywhere. Especially it should not include black strings off the brane. Naturally, the e solution is valid in the region 
r <C ^. It includes the region close to the horizon which is essential for horizon studies. However, it is not valid in 
the asymptotic region r — > cx) and one is unable to identify the mass or satisfy the junction conditions on the second 
brane. In order to overcome these subtleties we use the linearized solution as an asymptotic boundary condition. 
Linearized gravity is valid in the region r ^ fj,. We are interested in the case fj, <^ £, therefore, there is an intermediate 
region fi <^ r <^ £ where the e solution and the linearized solution coincide. In this intermediate range we take the 
large distance behavior of the e solution and compare it to the short distance behavior of the linearized solution. As a 
result, the long range characteristics (mass identification and second brane conditions) arc introduced into our short 
ranged e solution. 

In section Hill we analyze the properties of the horizon. Clearly, the MPS solution, (zero order in e), includes a 
regular horizon with a well defined temperature. The horizon is modified in the e solution. Although the e solution 
is valid in the region close to the horizon, it is not obvious that the new horizon is regular. The e solution acquires 
a null surface defined by gu — 0. The normal vector is null and the static Killing vector is null, i.e. the surface is 
a Killing Horizon. The surface gravity of the horizon is not a constant which indicates that the horizon is singular 
^j. This leaves us with the following possibilities; (1) The horizon is singular. Which means that there are no small 
black holes in Randall-Sundrum scenario. (2) The e expansion does not converge on the horizon. It is valid only at a 
distance e away from the horizon, and therefore one cannot use it to calculate the surface gravity. (3) The boundary 
conditions should be changed thus changing the entire solution. The issue of the regularity of the horizon is under 
study and the detailed discussion is left for a future publication. 

II. AN EXPANSION CLOSE TO THE HORIZON 

The goal of this paper is to find a method to explore the horizon of a small black hole in RSI scenario. The 
configuration is characterized by two length scales; the five dimensional Schwarzschild radius fi ITeV^^), and the 
AdS curvature length £ (~ lOTeV""'^) 15]. The only dimensionless parameter is the ratio between the two e = ii/£. 
For a small black hole e is smaller than 1 and the metric can be expanded in e. 

The zero order metric is just the MPS solution 

ds^ = -(^1-1^^ dP + (^-=r-^ <^p'' + /'"'^V'' + sin^ ^dnl . (2) 
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Where we use spherical coordinates in the (4+ l)-dimensional bulk (t, p, tp, 9, (f>). For convenience we rescale the radial 
coordinate p and the time coordinate t hy p thus making them dimensionless 



p t 2- 

P = - ; * = - ; 9AB = p 9AB 
p p 



(3) 



The TeV brane is located at ■0 — 7''/2 and the direction perpendicular to the brane is ip — Q. The Planck brane is 
located at pcos?/' ~ It is outside of the region where the e expansion is valid. In general, after rescaling the 
coordinates by p, the ansatz for the metric is 



(1 — epcoaipY 



-2 ^ Vn{p, ij)e''dpdi; + ( 1 + 5] Un{p, J + sin^ i^d^ 

n=l V n=l / 



(4) 



The conformal factor in front of the metric is just the Randall-Sundrum conformal factor and was introduced for 
convenience. The rescaled equations of motion include the Einstein equations [2^ 



1 2 
Rab - ■^R9AB = 6e gAB 

and the rescaled Israel junction conditions [^J] on the TeV brane 



(5) 



(6) 



The Planck brane is located at a distance ~ £ from the point mass and therefore the junction condition on the Planck 
brane can not be evaluated. 

Clearly, zeroth order in e corresponds to vacuum Einstein equations with no branes. The zeroth order solution is 
therefore the MPS black hole 0I- First order in e turns the branes on but the cosmological constant of the bulk is 
still absent. This case corresponds to an empty bulk with a brane and a point like mass on the brane, it is the main 
subject of this paper. Second order in e turns on the bulk cosmological constant. 

A. First Order in e 

The bulk equations © in first order in e are still the vacuum Einstein equations. However, the solution will be 
different from the MPS solution since the brane tension is of order e and it will effect the solution through the junction 
condition ©. 

The ansatz for the metric, up to first order in e, is 



(1 — epcos?/')^ 



eBi dt 



p^ - 1 



(Ax dp 



+ 2eVidpdij + p2 (1 + eC/i) d^P^ + p^ sin^ ipdnj 



(7) 



We choose to work in a gauge where t, 9, 4> are kept unchanged, and g, 
equations is 

_ 2p^ - 1 , P^cot^ 

The gauge ^ can be achieved using a coordinates transformation of order e 

p — > p{l — eF) ; —>■(/; + eFtan-f/; 



p^ sin ^. A convenient gauge in the bulk 



(8) 



(9) 



The bulk equations are indeed invariant under the transformation 0. However, the junction conditions are not 
invariant and the function F{p, ijj) is not completely arbitrary. It is subject to junction conditions on the brane 



(10) 
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Using the gauge the bulk equations can be solved in terms of two functions. The gauge function F(p, and the 
wave function H{p,'ip). The solution is 



Bi 



8p cos-0 



4p2(3i/(p,^)+tanV'ff^(p,^)) , Up^p^ - 2)Ho{p) , 2p(2p^ - 



ip^-1) 



(p2 - 1)2 



(p2 - 1) 



p2 - 1 



(11a) 



(lib) 



[6p [h{p, ij) - Ho{p) ] - (2p2 _ 1) ^Hjp, i,) - %[p) ] - pF(p, ^) - pcot ^F^(p, t/;)} , (11c) 



2 tan^ -0 

~2p2(2p2 - 1) sinV' - (4p2 - 1) isaii>H'^{p) + p(2p2 - 1) iimi>H'^{p) + 2^^ tan?AiJ,p(p, i/') 

p'F,^{p,4>) 



^ pi2p^ ~ 1) tan^ 0i?,^(p, V) 
where i?o(p) = H{p, ^ 



p tani/'i^^p(p, V) 



— 1 p^ — 1 

The wave function H^p^ip) is subject to the differential equation 



(lid) 



(p^ - 1) ( H,pp - ii/.p 



H 



sm ?/' cos ■(/' 



(12) 



The general solution can be casted as a combination of Associated Legendre functions in p and Hypergeometric 
functions in tp 



Hip,i^) 



dXR{p; A)1'(V'; A) 



i?(p;A) = pv/p2 - 1 a(A)gi^_i)/2(2p2-l) + &(A)P(\_i)/2(2p2-l) 
*(^;A) = c(A) 2Fi( 



1-A 1 + A 5.2m diX) ^,-2-A -2 + A 1.2m 

-^,^^,7:,sm ^2i^i( — ^ — , ^ — ,-7;,sm V) 

2 z z! sm tp 2 2 2 



The first hypergeometric function in Ea. (|15|l can be expanded by elementary functions as 



,1-A 1 + A 5 . 
"'^~'2' 



2^i(^^,^^,^,sin2 0) 



3[A cos(AV') sin(2V') - 2 cos(2V') sin(AV')] 

2A(4- A2)sin^ V 
3[sin(27/;) - 2V;cos(2V;)] 

8 sin^ ip 
3[4V'- sin(4V')] 
32 sin^ ip 



(13) 
(14) 
(15) 

(16a) 
(16b) 
(16c) 



The function H{p,tlj) is subject to the following boundary conditions. 

No black strings. The only source is the (point) mass on the brane at "0 = tt/2. Therefore, the solution must be 
regular everywhere. In particular one should eliminate the black string which lies in the bulk on a line perpendicular 
to the brane at = 0. The place to look for a string like singularity is the Kretchmann scalar (RabcdR^^^^)- The 
Kretchmann scalar depends only on the combination 3iJ + ta'mjjH_^. Using the expansion ()15|l one finds that 



S^i^P; A) + tan^*,^(0; A) - '^^^^^^^ + 0{i^°) . 



A solution free of black strings requires 

d{\) - d2(5(A - 2) 
The hypergeometric function associated with c?2 is simply 



2-^1 



-2-A -2 + A 1.2, 



= 1 



(17) 



(18) 



(19) 



A=2 
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Compatibility with Linearized Gravity. At large distances, p oo, the solution is given by Linearized Gravity which 
is an expansion in the mass, M ~ to first order. Iterating the linearized solution to include higher orders of M 
one finds that only integer powers of M (even powers of fi) enter the metric. Recalling that e = fi/£ ^ Af^/^ and 
p — p/ p, ^ M~^/^, it is clear that ep~" ^ (n+i)/2_ Since we allow only for integer powers of M, the functions H{p, ip) 
and F{p,ip) must be odd functions of p. The function R{p; X) [Ea. (|14|l ] has the same parity as p^^"^. Therefore, the 
integral over A must turn into a sum over even values of A. 

Junction conditions on the hrane. To first order in e, Israel junction conditions on the TeV brane are simply 



Si,v(P:f)=0; Ai,^(p,|)=0; V,{p,^) = Q 



Using the function these conditions imply that 



ttA 2 

dAi?(p;A)c(A)cos( — ) = --(p3_p) 



(20) 



(21) 



The asymptotic behavior of the Legendre functions of the first kind is: p\/ p^ — lP^\_^^y2(2p^ — 1) '--^ p^+-^. Therefore, 
to prevent higher powers of p, one must eliminate the Legendre functions of the first kind except for A = 2 which is 
the only function that can generate an asymptotic behavior p^. In particular 



&(A) =&2^(A-2) 



(22) 



Still, the Legendre function of the first kind, with A = 2, includes a term (log p)/p in the asymptotic expansion. In 
order to eliminate the logarithmic term from the condition (|21ll one must use the combination 



\=2 



This combination can be obtained by setting 

a(A) = 



n=0 



a2„<^(A-2n)-— 62'5'(A-2). 



(23) 



(24) 



Asymptotic boundary conditions. One would like to impose boundary conditions at p — )■ oo. The problem is that 
the e expansion is valid only for 1 < P ^ 1/e (which corresponds to /z < p <C ^). However, the asymptotic form of 
the metric is given by the linearized solution. The later is valid for 1 ^ p (/j, ^ p). As long as e <C 1, there is an 
intermediate region 1 ^ p ^ 1/e where both the e solution and the linearized solution are valid. Therefore, the large 
p limit of the e solution should coincide with first order in e of the linearized solution. The appropriate asymptotic 
boundary condition is 



lim pigAB-dAs) = • 

p — 'OO 

To summarize, collecting Eas. (tT^ . ((THll . (j^ . and (pi|l . the function H{p,4>) is 

3[sin(2i/') - 21/^ cos(2V')] 



(25) 



H{p,i;) = p^p^ - 1 



aog-i/2(V-l)- 



8 sin"* V 



?V2(2p'-l) 



02- 



3[4Vj - sin(4?/>)] 




32 sin^ V sin^ i/" 

2 5Q^_i/2(2p2 - 1 



3[8V' + 4V' cos(4V') 



dn 



3[4V'- 



327r2 sin^ V 
sin(4-0)] 



3sin(4V')]^ 



32 sin'' ip 



E l 2 , 3[ncos(2rt-0) sin(2V') — cos(2'(/;) sin(2nV')] 
«2«Q„-1/2(2P -1) «„n _.2Ae;.3 



8n(l - n^) sin"* ip 



(26) 



It includes a discrete set of parameters; ^2, ^2, and a2n {n ~ 0,1,..). These parameters should be fixed by the 
asymptotic boundary condition H25|l and the junction condition H21f) . 



6 



So, we turn first to the asymptotic behavior of the solution. The linearized solution appears in appendix 1X1 
Expansion of the linearized solution to first order in e appears in Eas. (|A23|l . The latter imply 



cos4> - 



sini/; 



5 3 V 

— - cos Tp — -i/jsinip + — — - — p cos 'ipWi „ 
4 2 Asm^p 

11 5 2-0 3 

— cos V' H — cos(3V-') H — sin tp + sin tpWi 

8 8 smyj 2 



cos2^/)(2- 5cos(2-0)) 3V'cos(3-0) 1 . 



1 



Ssin^ 



. 2 , + TTPsinV'VFi^p - -cosV'VFi,^ 
sm ip Z Z 



(27a) 
(27b) 
(27c) 
(27d) 



The function Wi{p,^) is arbitrary, but is subject to the boundary conditions 

W^i(p,|)=0 ; W,Mp,^)^^. (28) 

To simplify the analysis [to avoid the function M^i(/9, -0)], one should look at the following combinations; Si, Ci = 
(Ai/(cos 0)) + (pC/i/(sinV')) „, £>! = sin^ V^i + cos^ ipUi + Vi sin(20)/p, and C/f I , , 



COS V' ■ 



1 



-I- 



i/'=7r/2 



psin^(20) 
5 

IGpsini/; 
. 



20cos(2'0) - sin(2V') + - sin(40) - 3-0 



sin(4'0) - 40] 



(29a) 

(29b) 

(29c) 
(29d) 



For the purpose of comparing the e solution and the linearized solution, the e solution should be expanded in inverse 
powers of p only up to order ep^^ ^ M. Expanding the solution (|lll I26|l to first order in p"^, and casting the gauge 
function as F{p, -0) = pFi{tp) + i^2(0)/p, the asymptotic behavior of the e solution is 



CI 

^lli/'=7r/2 



3^ 
4 
1 



psin^(20) 



00-2^^1(0^) 

sin i^Fi (V-) + 8 cos sin^ il;F[ (0;) 



962 , . 



-Snao (30cos(20) — 3 sin 0; cos + tt sin^ -0) + ^— (sin(4?/') — 40;) 



—p cos sin^ ip[n — 362] H — 

TT p 

362 

' Stt 



(0^ — 20 cos(20) + sin cos — tt sin'^ 0) 



COS0 + 5cos(30) 



37r^ao 



+ Fi'{n/2) 



60 
sin0 
1 
P 



4 sin0 

+ 2sin2 i/'Fi(0) + sin(20)F{(0) 



STT^ao 277r2a2 97rd2 



32 



64 



4562 
16 



+ F.^'{7r/2) 



Comparing and Bf [Eas. H29a|l . (|30a|l ] and using the junction condition for F{p, 0) [Ea. pU|) ] one finds that 



COS0 ■ 



TT 

sin0 2 



Going next to and L'f [Eqs. if^ . (HE))]. One finds that 

O2 = 



Finally, Eas. pMj) and pMjl imply that 



J^2"(|) = ^(56 + 144^2 + 27^02) . 



(30a) 

(30b) 

(30c) 
(30d) 

(31) 

(32) 
(33) 
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The rest of the set a2„ can be evaluated by means of the junction condition on the brane (|21(l 

a2„(-l)>VP"-iy;-i/2(^P'-l)- ^ 

n=l 



J2 a2„(-l)>V7^Qi-i/2(V - 1) - — Pv/7^Q-i/2(2p' - 1) 
3/'V^(^i^/2(2p^-l)+;^ k ■ 



(34) 



To simphfy Ea. (|34|l . we set z 

E«2„(-im\-i/2W 



— 1 and divide the equation by — 1 

2 



(35) 



according to a theorem by Banerjee the toroidal functions Ql^_i^2(^) form a complete set on (1, oo). The function 
on the right hand side goes to zero as p ^ oo. Therefore it can be e xpan ded in a series of toroidal functions of the 
second kind i.e. the left hand side. The prescription given by Banerjee |23| is quite complicated and we will not follow 
it here. Alternatively, the coefficients can be found recursively from the asymptotic expansion of Ea. (|34|l 



02 



"3^ = " 9^ ' 



41 
2br 



(36) 



A single parameter d2 remains undetermined. 



III. THE BLACK HOLE HORIZON 



The existence of an event horizon is a delicate issue. It is clear that before turning on the brane (e — 0), there 
exists a single horizon at pn = 1. There is no guarantee that first order in e will not ruin the horizon. 

We assume that the new horizon is located at p^f = 1 + eC- We expand gtt{p = 1 + eC) = to first order in e to get 



2eC + 2e 



= 



(37) 



with f3 = (6(^2 + l)/8 depend on the parameter c?2- 

The gauge function F(l,-!/)) introduces a ip dependence. But, there is still a gauge freedom to choose F{l,ip) — 
which is in complete agreement with the junction conditions (I10|l and (|31I33I) . We choose the gauge F{1, t/j) — from 
now on. 

The logarithmic term comes from the associated Legendre functions of the second kind which close to 1 behave like 



py7^QLi/2(2p' -1)^-2 + ^ [- log |P - 1| + Cn] + 0{{p 1)2) . (38) 

This is the first indication that the solution (|26|l is not regular on the horizon. The second indication comes from the 
surface gravity. The surface gravity calculated from the solution 12()|) is not constant on the horizon. According to 
fl9j this shows that the horizon is singular. 

There are few possible explanations. One option is that the configuration indeed becomes singular already at first 
order in e. This means that there are no small black holes in Randall-Sundrum scenario. Another possibility is that 
although the e solution is valid close to the horizon it does not converge on the horizon. This can be seen if one tries 
to expand g^^^ in e, clearly, the expansion does not converge on the horizon {gu = 0). 

The regularity of the horizon is still under study. The detailed discussion of this issue is left for a future publication. 



IV. SUMMARY 



Randall-Sundrum 1 model makes the gravitational interaction strong on small scales and might allow for small black 
holes on the TeV brane. Interesting phenomena might arise in the presence of TeV scale black holes. Especially, the way 
such primordial black-holes accrete or radiate matter when immersed in a hot plasma would affect the cosmological 
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evolution of the universe. The behavior of a growing or evaporating black hole is related to its thermodynamic 
properties. Thus the study of its horizon is required. 

The purpose of this paper is to introduce an approximation method for studying the horizon of a small black hole 
in the Randall-Sundrum I scenario. 'Small' black hole means that the five dimensional Schwarzschild radius of the 
black hole is smaller than the curvature length of the bulk AdS manifold. The method suggested in this paper is to 
expand the metric in the dimensionless parameter e = \/G^M ji, and to fix the asymptotic behavior using the weak 
field approximation (linearized gravity). Zeroth order in e is just a Myers-Perry flat background black hole. First 
order in e corresponds to the presence of a brane with brane tension which is of order e in a flat bulk. Second order 
in e includes also the bulk cosmological constant which is of order e^. 

We fix the asymptotic behavior in the following way. The e solution is valid in the region r <ti £. It includes the 
region close to the horizon which is essential for horizon studies. However, it is not valid in the asymptotic region 
r — > oo and one is unable to identify the mass or satisfy the junction conditions on the second brane (at distance ~ £). 
In order to overcome these subtleties we use the linearized solution as an asymptotic boundary condition. Linearized 
gravity is valid in the region r ^ /i. We are interested in the case /i <C ^, therefore, there is an intermediate region 
fj. <^ r <^ £ where the e solution and the linearized solution can be matched. In this intermediate range we take the 
large distance behavior of the e solution and match it to the short distance behavior of the linearized solution. As 
a result, the long range characteristics (mass identification and second brane conditions) are introduced into the e 
solution. 

We should mention that the method of expanding in e and comparing with the linearized solution as an asymptotic 
condition is applicable to other scenarios as well. For example: Randall-Sundrum single brane scenario, de Sitter 
bulk (expanding in the ratio between the Schwarzschild radius and the curvature length), compact extra dimension 
(expanding in the ratio between the Schwarzschild radius and the compactification radius), etc. 

In this paper we have found a solution up to first order in e. We require that the curvature singularity resides only 
on the brane. No black strings are associated with our solution. The mass of the object consists of matter completely 
confined to the brane. 

Although the e solution is valid in the region close to the horizon, it is not obvious that the new horizon is regular. 
The e solution acquires a null surface defined by gu = 0. The normal vector is null and the static Killing vector is 
null i.e. the surface is a Kil ling Horizon. But, the surface gravity of the horizon is not a constant which indicates 
that the horizon is singular [T3. This leaves us with the following possibilities; (1) The horizon is singular. Which 
means that there are no small black holes in Randall-Sundrum scenario. (2) The e expansion does not converge on 
the horizon. It is valid only at a distance e away from the horizon, and therefore one cannot use it to calculate the 
surface gravity. (3) The boundary conditions should be changed thus changing the entire solution. 

The issue of the regularity of the horizon is under study and the detailed discussion is left for a future publication. 
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APPENDIX A: LINEARIZED GRAVITY 

1. Randall-Sundrum I Background 

RSI scenario consists of two flat 3-branes embedded in five dimensional AdS bulk. The equations of motion in the 
bulk are the five dimensional Einstein equations 20] 

1 6 
Rab - -^RgAB = SttGsTab = -pdAB ■ (Al) 

With an additional Z2 symmetry with respect to the brane, the Israel junction conditions are 

±6 

2 {K-f^, - K^,) = 8ttG5S^, = —7^, . (A2) 

The upper (lower) sign corresponds to the negative (positive) tension brane. 

Parameterizing the bulk with the coordinates = (a;^, w), where w is the coordinate perpendicular to the brane, 
a conformally flat ansatz for the metric is 

dsl = a'{w)j^ABdy^dy^ . (A3) 



9 



We choose a coordinate system in which the TeV brane (with negative tension) is located at w = 0, while the Planck 
brane (positive tension) is at — l(\ — A). A is a very small number. It is the ratio between the TeV scale and the 

Planck scale, namely A ^ 10"^^. The conformal factor in ljA3|l is c?{w) — ^ ^_^^| ^ . It has the value 1 on the TeV 

brane (at w = 0) and the value A^^ on the Planck brane. In this scenario the standard model particles are assumed 
to be confined by some mechanism to the TeV brane. 

Let us assume now that there is a matter distribution on the brane such that it is confined to a small region on 
the brane at w = 0, therefore, the energy momentum tensor is non vanishing only at (r < rg, w = 0). However, if we 
are interested in the solution far from this region, we can use the linearized version of Eas. ljAllA2|l . The perturbed 
metric is given by 

d&\ = (w) \^AB + hAB {x" , w)] dy^dy^ . ( A4) 

The perturbation Hab is assumed to be small. The source for the perturbation is the matter on the brane, characterized 
by the energy momentum tensor, T^i^. The total energy momentum tensor on the brane (which enters Israel junction 
condition iP^ l is 5^1. = a^fj,^ + T^i,. 



2. Conserved Momentum 



The definition of a conserved momentum is usually done in an asymptotically Minkowski n + I dimensional space- 
time as = J d^^xT"^^ . It is associated with energy momentum conservation T^g — 0, where the use of ordinary 
derivative instead of covariant derivative is essential. 

In RS models, one deals with asymptotically conformal Minkowski space-time, and the energy-momentum conser- 
vation condition is T"^^ -b = 0. One can verify, however, that the energy momentum conservation condition can be 



written as 



Tl^^ = 1 (yr^Tf ),5 - ^^tE5a^ = . (A5) 

Therefore, since the conformal factor a{w) depends only on the fifth coordinate, one is able to define a conserved four 
momentum P^^' = J d^xdw ^—gT*^. One can verify that the momentum is conserved using 



p(5) 



j d^x dw [(T^^T; ),, + iV^rp,^] = - \im^ J dw dn r^n.^T^ - J d^x dw (V^TJ^'),^ . (A6) 



The first integral on the right hand side vanishes provided the matter is confined to a small region. The second integral 
on the right hand side is treated differently in different scenarios; in RSII, with a single brane, the w integration is 
over the interval (—00,00) and vanishes provided the matter is confined to the brane and the energy-momentum 
component T™ is continuous. In RSI with two branes the only requirement is the continuity of TJf . 

It is worth noting that if the energy-momentum tensor is traceless (for example four dimensional Maxwell field), 
then one is able to define a conserved five momentum P^^-* — J d^x dw ^J—gTj^. In other words, the presence of the 
trace brakes a symmetry and therefore the theory has only a conserved four momentum. 

The solution to Eq.ljXS)) is given by yf^T^ = doQl^ , with Q^^ = -Q^^- The exact form of Q^^ can be 
calculated from the linearized Einstein equations. The conserved momentum can be calculated using 

P^5) = j d^xdwd^Q'^° + J dndwr'^n.Qf . (A7) 

The first integral has contributions only from the discontinuities of Q^^. This can be evaluated using the linearized 
junction conditions and turns out to be exactly the four momentum of matter confined to the branes 

/ d'-dwd^Q-' ^ I d^x rO|br,,es = ^i'^ ■ (A8) 
The five dimensional momentum can be summarized as 
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The energy Pq is given by 

Pp(5) ^ pii) ^lij^ _i_ y a3(u;)r2n' [-/io^,o - /^l, + + h^^,,] . (AlO) 

If the only sources for the perturbation reside on the branes, then the five dimensional and the four dimensional 
momenta are equal, and the integrals in Eas. (jA9IA10|l must vanish. 



Linearized Solution 



Choosing the gauge ggg = a^{w)r'^, the metric perturbation caused by a single point-like mass, A/, located on the 
TeV brane is given in radial coordinates as 



htt{r,w) = 



2C 2G5M 



hrrir, w) 
hgg{r, w) 

hww{r, w) 



r 3r£ 

2C 2G5M r ~ 

j 2_ 

r 3r£ 



2Fi(a,/3) 



W{r, w) 



Fi{a,P)-aFi,a 



= , 



= -2 



fW{r, 



V i-w 



3£2 



i — w \ £ — w 



(Alia) 

(Allb) 
(Allc) 
(AUd) 

(Alle) 



The term — is the contribution of the zero-mode. The constant C cannot be determined by the bulk equations nor 
by the junction conditions. The only way to determine C would be through the definition of the total mass of the 
configuration given by Eq.JUni. 

The terms with Fi , F2 are the massive tachyonic modes of the linearized equations. These functions depend on the 
dimensionless variables a = r/£, /3 = {£ — w)/£, and are given by 



dz sin(zQ;) 



Ki{z\)h{z) - h{z\)Ki{z) ' 



F2(a,/3) dz 



sin(zQ;) 



sin(zQ;) cos(2;a) 



K,{zX)lAzP) - h{zX)K,{zP) 
Ki{zX)h{z) - hizX)Ki{z) 



(A12a) 
(A12b) 



W{r,w) is a gauge function which cannot be determined by the bulk equations. It represents a coordinates transfor- 
mation of the form 



r —f r ; w w + W{r, w) 
It is subject to the following boundary conditions 

G5M 



(A13) 



W{r,w = 0) = 



; VK,„(r,w = 0) = - 



G5M 



3r ' ' 3r£ 

W(r,w = ^(l-A)) =0 ; W,Ar,w ^ £{1 - X)) = 



(AM) 



On the TeV brane {w = 0) W{r,0) behaves like a massless four dimensional scalar mode, and is called "the radion". 

The definition of the mass and the hierarchy should be seen at the asymptotic region r ^ £. The asymptotic 
expansion of the functions Fi , F2 is given by 



Fi(a» 1,13) 
i^2(a» l,/3) 



2A2 



1-A2 ' 

W - A^) 
a2(l- A2) 



(Al5a) 
(A15b) 



Using Eg. ljAlOjl and assuming that there are no other sources besides the point mass on the TeV brane, one can fix 
the constant G in the zero mode G = The zero mode tends to be repulsive. 



11 



The gauge invariant force acting on a slowly moving test particle is given by 

2 t — W 



(A16) 



The gravitational potential along the TeV brane at distances much larger than the curvature length r ^ ^ is given 

by 



V = - 



2GsM\^ G5M 



3ri{l - A2) 



3H 



(A17) 



The first term is the contribution of the zero-mode and the massive modes. These sum up together to give the correct 
hierarchy with G4 = 3£(j^i'\2) • 

The last term in Eq. (|A17p is the contribution of the radion and is clearly destroying the hierarchy and producing a 
long range strong (TeV) gravitational attraction. In order to save the hierarchy it is necessary to include a stabilization 
mechanism that removes the radion at large distances. A stabilization mechanism was suggested in (2Sj. It introduces 
a bulk scalar field which is coupled to the radion by some potential and cause the radion to acquire mass. Thus, the 
radion acquires a gravitational potential of the form e~'"''''/r and does not contribute at distances r ^ ^ log(-p-). 
In general, the introduction of a stabilizing bulk scalar field will change the RSI background, namely the conformal 
factor, and might have an impact on the graviton modes as well. 

In this paper we are interested in the short range region fj, < r < £. Therefore, assuming that the radion becomes 
massive by some mechanism but is still light compared to the curvature length < j 10— lOOGeT^, we can assume 
that in the region r < £ the above discussion is valid and the metric for the linearized solution is given by Eas. (|IlT|l . 
In addition, we assume that the stabilizing bulk scalar field is massive enough not to interfere in the region /i < r < £, 
thus assuming 

The asymptotic expansion of the functions Fi , F2 in the region r i is given by 



f'2(«< l,/3) 



7r(a2 + (1 - /3)2) 
4a(l - /3)/33/2 



^(^-^)^arctan( " 



(3/32 + 2/3 + 

7r(a2 + (1 _ /3)2)2 + 4^Q, 

Using these expansions, and also expanding the gauge function as 



1-/3 



) , 



f3 



1 , a 

— arctanf - 

a 1 — P 



(Al8a) 
(A18b) 



W{r,w) = ^ 



Wo{r, w) + jWi{r,w) 



one can write the expanded form of the metric IjAlip , to first order in 1 / 



htt{r, w) 

hrr{r, w) 
hgg{r,w) 

h^^{r, w) 

hrn,{r,w) 



2 7 



(r2 + tiP'Y 



3 ,r 3u) 

— arctan(-) - — ^- — ; 

3 , r 3u'(w2 

-— arctan(— ) + 
4r w 4(r"^ -\ 



+ Wo{r,w) 
-r2) 



+ rWo^r{r, w) 



= 



-M^Wo,tt,(r, w) 



1 



Wo(r,w)+Wi.^{r, w) 



w(3 



r 



■ arctan( — ) + , , , 
4r2 w 4r(»-2 



w 



1 



8^ 



arctan( — ) 



2')2 

3w2(3r2 + w2) r 



8r(r2 + ^2)2 



The boundary conditions for W{r, w) ljA14p are 



(A19) 



(A20a) 

(A20b) 
(A20c) 
(A20d) 



(A20e) 



M^o(r,0) = - 



Wi{r,0) = 



n 
4r 



(A21) 
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In order to compare the linearized solution (|A20|I and the e expansion dJ, one has to perform a coordinate trans- 
formation to rescale the linearized-mctric (IA20|I 



I /''' \ 1 M 

ii = arctan( — ) ; qab -^gAB e = - 



After transformation IIA22II the linearized metric (IA20II is 



(1 


- epcosV')^gft 


(1- 


- ep cos iP)^g^p 


(1- 


ep cos ip)'^g^^ 


(1- 


-epcosiPfg'^^, 


(1- 


-epcosipfg^g 



1 e 

+ — -- 
P P 

1 e 

— + - 

P P 



cos i) - 



sin-f/) 



5 3 "0 

- - cos ip — -ip sin -0 + "T^ — ■< — p cos iJjWi „ 
4 2 Asimp 



11 5 21/) 3 

— cos -0 H — cosfSV') H — -0 sin ■(/; + sin tpWi ^ 

8 8 smi/) 2 

cos2 0(2 - 5cos(20)) 30cos(3V') 1 . 



Ssin-!/; 



n • 2 , + TTpsinV^VFi^p - -cosV^W^i.^ 
8 sm tp 2 Z 



p^ siv? -0 . 

The first term of the gauge function W{p, tp) (|A19p is fixed by setting e = 



2fj.p 



sintp 



(A22) 

(A23a) 

(A23b) 

(A23c) 

(A23d) 
(A23e) 



(A24) 



The second part of the function W{p, ip) IIA19|) is yet undetermined, but is subject to the boundary conditions 



(A25) 
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